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Abstract 

We compute the structure function g2 for a gluon target in perturbative 
QCD at order as- We show that its first moment vanishes, as predicted by the 
Burkhardt-Cottingham sum rule. 
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A computation of the structure function g2 for a target quark at order as in 
perturbative QCD has been performed in ref. The interest in that computation 
was mainly driven by the possibihty of performing a direct test of the vahdity of the 
Burkhardt-Cottingham (BC) sum rule in perturbative QCD. The BC sum rule 
states that ^ 

/ dxg2{x,Q^) = {) (1) 

on the basis of general arguments about the analytic structure in the complex v = 
Q^/(2x) plane of the amplitude whose imaginary part gives g2 (a detailed discussion 
on the derivation and the validity of the BC sum rule is given in ref. [§]). It is therefore 
interesting to check whether eq. (|I|) is valid in perturbative QCD. The result of ref. 
(later confirmed in ref. Q]) is that the first moment of g2 actually vanishes at order 
as for a massive target quark (at leading order, g2 itself vanishes for a target quark). 

In this note we present an analogous calculation for a gluon target. We will show 
that also in this case the BC moment vanishes in a wide class of regularization schemes 
for the coUinear singularities. 

Structure functions are defined starting from the Fourier transform of the forward 
matrix element of the product of two electromagnetic currents between polarized 
states: 

W'^''{p,q,s) = ^J d'xe"^^{p,s\J'^{x)rmp,s), (2) 

where p and s are the target momentum and spin four-vectors, respectively. In the 
case of a target gluon, the tensor is given by 

W'^''{p,q,s) = W>^''P'^{p,q)e^el, (3) 

where the gluon polarization vector e carries the dependence on s. Only the antisym- 
metric part of W^'^ is relevant for the computation of polarized structure functions; 
we have 

where in the last step we used the symmetry of W^^'^" under the simultaneous ex- 
changes fi ^ u, p ^ a. We are therefore interested in computing the antisymmetric 
part of the gluon polarization density matrix pp^. = CpC*. To do this, we assume that 
the gluon in the initial state is off the mass shell, p"^ ^ 0; this allows us to define a 
longitudinal spin vector for the gluon: 

= AA^ I - — g° 1 , (5) 
V PQ J 
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which satisfies the transversity condition = 0; A is the gluon hehcity, and the 
normahzation factor (real and positive) is related to by 

o2 



where 



P 



p2q2 



'1 + 



(6) 



(7) 



with the usual definitions = — g^, x = Q'^/{2pq). We see that we can choose 
so that = 1, provided and have opposite signs. We can now write the 
antisymmetric component of e^e* clS cl function of p and s: 



the normalization is fixed by the condition that, for ^ 0^, the imaginary part of 
Pi2 is equal to —A/2, where A is the gluon helicity. 

The structure functions gi and g2 are conventionally defined by means of the 
following, general parametrization of iW'^: 



iJ\p'^ 
pq 



qs 



gi{x, Q )sa + g2{x, Q ) [sa Pa 



pq 



(9) 



where s is now a generic spin vector (not necessarily longitudinal, since the gluon is 
off the mass shell). One can define projectors P 



pV> 



1, 2 such that 



One possible choice is the following: 



1,2) 



PL = P-%..,q- (p'+'^-J-s^ 



s"^ pq 



with 



P = 2J\p'^\qs 



p2q2 
{pqf 



{qsf 
q2g2 



(10) 



:iii 



(12) 



(13) 



Notice that P = even for 7^ if s is purely longitudinal, as one can see using 
eqs. (|^) and (^); we must assume that s has a transverse component until we have 
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projected out the structure functions, which are by construction independent of s. 
Collecting everything together, we finally obtain 

g,{x, Q') = -p-h^'^^' (^ps + ^^ss^ -4^= e^'-^^p^spW.^.^ip, q) (14) 
g2{x, Q2) = _p-ie/^-7<5p^ (^q^ _ ^ss^ e^'-'^^p^SpW^^.^ip, q). (15) 



Even before explicitly computing W^upa^ we can check that gi and g2 given by eqs. (0) 
and ( p!5|) are indeed independent of s, as they should. In fact, it is easy to prove that 
only a term proportional to gfipP^qa gives a non-zero contribution to gi, while in the 
case of g2 the only surviving term is proportional to g^pq^qa- Inserting these terms in 
eqs. (0) and ([TsD , the s-dependence is seen to cancel against the factor P^^. 

We now proceed to compute W^yp^ip^q). The calculation in this case is much 
simpler than in the quark case, because at order as there are no loop diagrams that 
contribute to the relevant amplitude. For this reason, no ultraviolet or soft divergences 
are involved. We have 

= i^]^ E / d<P^'^ (A« + 45) (A« + Afy)\ (16) 

c colour 

where A^c = 3 is the number of colours, and is the two-body phase space. The 
amplitudes A^^\ A^'^'^ correspond to the diagrams of fig. The singularities that arise 
when an on-shell gluon radiates massless quarks in the coUinear configuration may 
be regularized either by a non-zero gluon virtuality p^, or by a non- vanishing quark 
mass m. We will keep both p^ and m different from zero at this level, and we will 
discuss later the behaviour of our results in the limit p^, m? — > 0. In the photon-gluon 
centre-of-mass frame, the momentum k of the produced quark is 

E 

k = —{l,0,(3sme,(3cose), (17) 

where 

E^ = {p + qf=p' + Q'^ (18) 

X 

and 



The two-body phase space dcp^'^^ takes the form 



dcpi^) = J-dcose. (20) 
IQtc 



Figure 1: Diagrams contributing to W^'^ for a gluon target. 



The denominators of the virtual quark propagators appearing in the amphtude are 
given by 

{k- qf -m^ = -pq (l+ (3 (3 cos e) (21) 

{k-pf -m^ = -pq {l- (3(3 cose) . (22) 

The phase-space integration is therefore singular for cos^ ^ ±1 when /? = /3 = 1, or 
equivalently p^ = m? = 0. The calculation is straightforward (we have performed it 
with the help of the algebraic manipulation program MACSYMA) ; the cos 6 integra- 
tion can easily be performed by observing that, after inserting eq. (|T^) in eqs. (|I^) 
and ([T5|), the numerator of the integrand expression is a degree-2 polynomial in the 
invariants kq, kp and ks. Terms proportional to powers of kq and kp can be integrated 
immediately, since their dependence on cos 6' is given explicitly by eqs. ( [21| ) and (|2^). 
Terms containing powers of ks can also be expressed in terms of integrals containing 
only kq and kp. Consider for example 

p = J d(P^^'^ f{kp, kq) P, (23) 

where f{kp, kq) is a generic scalar function. The result must be a linear combination 
of q'^ and p'^: 

It" = Aq^" + Bpf", (24) 
and the scalar coefficients A and B can be obtained by solving the system 

Aq^ + Bpq = q^I^ (25) 
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Apq + Bp"^ =Pf,I'', 



(26) 



so that, finally, 



s^F = J dcj)^"^^ f{kp, kq) ks = Aqs. 



(27) 



Terms proportional to (A;s)^ can be treated in a similar way. Therefore, all phase-space 
integrals are of the type 



dcose{i + pp cos ey{i- pp cos ( 



with a and b integers between —2 and 2. We obtain the following results: 



(28) 



9i 



e^ac 1 



P 



^ {Ap^x' - Sp'P'x' - 8p'x' + 12x' - 2p'P^x 



-6p^x + Sp^p^x + 12p^x - 12x -p^ + 2p^p^ + 3/3^ - 2p^p^ - 5/3^ + 3) 



-^{iP^x^ - SP^x^ + 12x^ - 4p^x + Up^x - 12x -(3^ + 3(3^ - 5(3^ + 3) 



(29) 



92 



e^as 1 
' 87r p^x 



P 



P^P^x^ + AP'^x' - \2x' + 2^'P''x 



1 - P^P^ 



+2p^x - SP'P'x - SP'x + 12x - 2P'P^ -P^ + 2p'p' + 4p' - 3) 



-^{4p^x^ - 12x2 - 8p^x + 12x - + 4/3^ - 3) 



(30) 



where e is the electric charge of the produced quark in units of the positron charge, 
as is the strong coupling, and 



log 



1+PP 
1-PP 



(31) 



The collinear singularities are collected in the factor L, which diverges logarithmically 
when both and p"^ go to zero. The structure function gi was first computed in 
ref. 1^ for m = 0, p^ < 0. The general case ^ 0,^^ ^ was considered in 
ref. 0. Our result for gi, eq. (pQ]), is different from the analogous formula obtained 
in ref. 0] . The origin of this discrepancy is the fact that the operator used in ref. 
to obtain gi from W^'^ actually projects out the desired structure function only in 
the limit p^ ^ 0, while is kept non-zero elsewhere at this stage. However, the final 
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result of ref. is correct in the physically interesting limit, as we shall see later. 
Equation (pO]), on the other hand, is a new result. 

For Q'^ ^ oo with 



fixed, we find 
91 

92 



Arx^ — Qrx^ + 2rx — 4x + 3 
rx^ — rx — 1 

6rx^ — lOrx^ + 4rx — 4a; + 3 



rx^ — rx — 1 



(2x - 1)L 
- {2x - 1)L 



In this limit, L takes the form 



Q 



L = log — - — log rx + 



X 



X 



(32) 

(33) 
(34) 

(35) 



It is interesting to notice that the terms proportional to L, which contains the coUinear 
divergence, cancel in the sum gi + g2- 

As a test of the correctness of our calculation, we can check that we reproduce 
the known results for gi. Indeed, eq. ( |33| ) gives 



0: gi 



47r 



-4x + 3 + (2x - 1) f log ^ - lo, 



X 



1 — X 



/ dxgi{x,Q^) = 0; 
Jo 



r ^ oo 



9i 



Air 



-Ax + 2 + (2x - 1) log - logx 



dx gi{x, 



An ' 



(36) 
(37) 
(38) 

(39) 



which coincide, for example, with the results of ref. . The fact that different choices 
of the regularization scheme lead to different results for gi (and in particular for its 
first moment) has important physical implications, and has been widely discussed in 
the literature M. 



We now turn to the structure function g2. From eq. (|34D we get 



0: g2 



r ^ oo : g2 



e^as 
Att 



4x — 3 — (2x — 1) log — - — log 



1 — X 



Ait 



6x - 4 - (2x - 1) log 



_p2 



logx^ 



(40) 
(41) 
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The first moment of g2 vanislies in botli cases; actually, it is easy to prove, using 
eq. (H), that 

f dxg2ix,Q^) = (42) 
Jo 

for all values of r. Therefore we conclude that the BC sum rule is satisfied also by 
the gluon contribution to g2 at order as, within the class of regularization schemes 
we have adopted. 

Having computed g2{x, Q"^), its n-th moment 

fi 



92 



/ dxx^ ^ g2(x,Q'^ 
Jo 



can be obtained for any n. In the two cases r = and r 
is given by 



(43) 



oo, the n moment of g2 



: 



oo 



92 



92 



in 



e^as 



n + l 
6 



3 

n 



1 



n 



4:71 



n + l 



4 

n 



n{n + 1 
4 



log^ - S{n) 



[n + l] 



2 

+ — - 



n 



nin + 1 



log 



_p2 



(44) 
,(45) 



where S{n) = J2k=i V^- Once again, we see that both expressions vanish when n = 1. 

The computation presented here is equivalent to the calculation of gluon coefficient 
functions in the light-cone operator product expansion of W^^'^; in that case, the 
quantities that are directly computed are odd moments of the coefficient functions; 
in the case of g2, only moments for n > 3 are obtained, and therefore no direct test 
of the BC sum rule can be performed. Such a calculation was performed in ref. |^ in 
the case m = 0, which corresponds to our eq. (^5)). The two results are in agreement. 

In conclusion, we have performed a calculation of the structure function g2 for 
a target gluon. We have considered various regularization schemes for the coUinear 
divergences, and we have found that in all of them the first moment of g2 vanishes as 
expected. 
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